Observing the evolution of a quantum system that does not evolve 
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This article deals with the problem of gathering information on the time evolution of a single 
metastable quantum system whose evolution is impeded by the quantum Zeno effect. It has been 
found it is in principle possible to obtain some information on the time evolution and, depending 
on the specific system, even to measure its average decay rate, even if the system does not undergo 
any evolution at all. 



I. INTRODUCTION 

A quantum system initially prepared in a state that is 
not an energy eigenstate will generally undergo an evo- 
lution. For very short evolution times, the transition 
probabilities towards other states increase as the square 
of the time, and a metastable quantum system, which 
is observed frequently enough, would remain essentially 
unchanged. This is the quantum Zeno effect, pictorically 
described as the fact that, in quantum mechanics, "a 
watched pot never boils". This effect has been subject 
to intense theoretical [J 0, H, 0, [H, 0] and experimental 
@, § investigations. 

The question this article wants to assess is whether 
any information about the time evolution, and, above 
all, the nondiagonal elements of the Hamiltonian that 
determine the transition rates, can be measured while 
the system is being forbidden to evolve by a series of 
projective measurements performed frequently enough to 
be in the quantum Zeno regime. 

At a first glance, the answer may seem to be negative, 
each measurement projecting the system wavefunction in 
its initial state and, upon measurement, all the informa- 
tion about time-evolved states of the system being lost. 

Moreover, if such procedure existed, it would give us 
the possibility to perform some quite paradoxical mea- 
surements when applied to normally irreversible dynam- 
ics. 

Given an unknown quantum system in an excited state, 
undergoing an irreversible spontaneous decay, we could 
be able to know if the system is stable or not before 
the decay takes place. We could even manage to gather 
information on its average decay rate, not only by mea- 
suring a single system instead of a statistical ensemble, 
but keeping this single system undecayed. 

We will show that, given the possibility to make the 
system interact coherently with a properly crafted envi- 
ronment, such measurement is indeed possible. 

The measured quantity will be a certain function of 
the coupling with all the possible final states; the actual 
amount of information we can obtain will depend largely 
on the specific system and how much we know about it. 

This work is organized in the following way: in section 
II, the general ideas exploited in the rest of the paper are 



introduced. These ideas will be used in sections III to lay 
down the full quantum mechanical calculations and the 
efficency assessment of the measuring procedure. Finally, 
in section IV, the developed theory will be applied to 
some specific model Hamiltonians. 



II. GENERAL FRAMEWORK 

Let us suppose we have a certain system, whose evo- 
lution is being slown down by a series of projective mea- 
surements. If the measurements are frequent enough, 
there is a high probability that each one of them will 
find the system in its initial state. In order to gather in- 
formation on the time evolution of the system, the only 
thing we can do, in the time interval between two succes- 
sive projective measurements, is to make it interact with 
some other measuring apparatus. 

The system, found in its initial state by one of the re- 
peated projective measurements, is left to evolve freely 
for a certain fraction of time before the following projec- 
tive measurement; then we make it interact coherently 
with some ancillary quantum system, such that the state 
of the system is recorded in the ancilla. Let us say that 
a will be recorded if the system is found in its initial 
state and a 1 if it is in some other state. Then the system 
is left to evolve again until the next projective measure- 
ment, that still finds it in the initial state. Now there are 
no correlations anymore between the system and the an- 
cilla, we can thus check the ancilla state independently 
from the state of the system. The state of the ancilla 
just after the projective measurement will not be a but 
a certain linear superposition of and 1, where the co- 
efficient of the 1 is generally not zero. It is due to the 
finite probability that the system, during the first part 
of the evolution (before the interaction with the ancilla) , 
evolves to some other state and, during the second part 
(after the interaction with the ancilla), evolves back to 
the initial state. 

Even if the projective measurements always find the 
system in its initial state, because of the quasi-reversible 
dynamics typical of the quantum Zeno regime, we could 
rely on these sort of higher order processes to gather in- 
formation about its time evolution. In order to give a 
more formal and general definition of a higher order pro- 
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cess, let us suppose that, during the time interval between 
two successive projective measurements, the system in- 
teracts several times with an environment, so that each 
interaction leaves in it a record of the state of the system. 
We will call a n-th order process a branch in which the 
system interacts at least n times with the environment 
and the environment records exactly n jumps, both from 
and to the initial state (fig. [IJ. 

According to this definition, in the previous example, 
the nonzero coefficient of the 1 is due to a second order 
process. 




FIG. 1: A schematization of higher order processes in the 
simple case of a two levels system. The system interacts four 
times with the environment at t = ti,t2,ts and t^. The pro- 
jective measurements are performed at t = and t = r and 
both of them find the system in its excited state. The solid 
lines show the trajectory of the system. The first one corre- 
sponds to a zeroth order process, with no jumps, the second 
and the third to a second order process, and the fourth to 
a fourth order process. Clearly, if in the end the system is 
to be found in its initial state, only even order processes are 
possible. 



Obviously, the probability of a n-th order process will 
generally be of the order of the n-th power of the prob- 
ability for the system to jump only one time, that is the 
probability of it being found by the projective measure- 
ment in a state that is not its initial state. This prob- 
ability vanishing in the quantum Zeno regime, it would 
seem impossible to gather any information on a single 
system while still blocking its evolution, ft will be shown 
that on the contrary, by using a carefully crafted interac- 
tion with the environment, instead of the simple one of 
the previous example, it is indeed possible to make the 
probability of the second order processes high enough to 
measure the desired quantity, while still finding the sys- 
tem in its initial state. 



III. FORMAL THEORY 

Let |e) be the initial state of our system and \gt) a 
set of mutually orthogonal states, orthogonal to |e), such 
that |e) and all the \gi) form an orthonormal base for the 
Hilbert space of the system. The index i will in general 
run over a discrete as well as continuous part. We will 



consider as ancilla a two levels system (a qubit) prepared 
in the state, for the moment arbitrary, \tp). Let H be 
the Hamiltonian of the unperturbed system and K(t) = 
e -iHt/n j£ g evo i u tion operator. The composite system 
is prepared in the state |y>(0)) = |e)|^>). A projective 
measurement is then performed on the system each r 
that, with a certain probability, projects the system state 
upon the state |e). We will suppose r to be small enough 
to allow us to neglect all terms of order greater than two 
in the time development of K(r). 

We let the system evolve freely for a time r/2, neglect- 
ing the auto Hamiltonian of the ancilla. The resulting 
state of the composite system will be 



\<p{T~/2)) = K(r/2)\eM 

= [K ee \e)+J2 K ie\9i)W, 

i 

where K ee = {e\K(r/2)\e), K ie = { 9i \K{t /2)\e) and 
K e i — (e\K(r /2)\gi) . Then we entangle the ancilla with 
the system by performing a controlled-M on it, condi- 
tionally on the system being or not in its initial state, 
that is we act on the ancilla with the identity operator 
if the system is found in its initial state and with the 
operator M otherwise. The operator M is for the mo- 
ment a to be determined unitary operator that acts only 
on the ancilla, whose two eigenvalues we will name e w+ 
and e^and the respective eigenvectors |+) and |— ). Let 
assume the time to perform the controlled-M to be much 
smaller than r, in order to be allowed to ignore the effect 
of the free evolution. We obtain 

\<p(r+/2)) = [K ee \e)+Mj2KMm- 

i 

After a second evolution of r/2, just before the next 
projective measurement, the state will be 



\ V (t-)) = K(t/2)Mt + /2)) 

= [(Kl + Mj2K ei K ie )\e) 

i 



The projective measurement will find the system in the 
state \e) with a probability 



Pe = \Ki 



I A Jiff 



-2Re[K 2 ee AK 2 e MMM, (1) 



where we use the completeness relation I = \e)(e\ + 
Ei \9i)(9i\ to write 

Y,K ei K ie = (e\K(T/2)(I-\e){e\)K(T/2)\e) 

i 

= (e\K(T/2f\e) - {e\K{r/2)\ef = AK 2 e . (2) 
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Supposing we have effectively found the system in its 
initial state, therefore the state of the composite system 
will be 



Taking \tp) as an equal weighted superposition of the 
two eigenvectors of M, that is \tp) = ^(|+) + we 
obtain from eq. [3] and [7] 



Wt + )) 



Kl + MAKl 



\eM = U eff \e)\iP). (3) 



The operator U e ff acts only on the ancilla's Hilbert 
space, so that we can say that, with a probability P e , 
the system is left unchanged and the ancilla undergoes 
an effective evolution U e ft. 

The probability P e is given by the expression 

P e = \Kl\ 2 +\AK 2 ee \ 2 +2Re[K 2 ee AK 2 e MM\M (4) 

where it is worthwhile noticing that P e , and conse- 
quently U e ff, depend on the initial state of the ancilla, 
due to the fact that the effective evolution operator comes 
from a non- unitary dynamics. 

Since we are interested in the Zeno regime, we can 
make a development of K(r/2) to the second order, ob- 
taining 



K(t/2) = 1 



itH t 2 H 2 



2h 



m 2 



T 3 H 3 

°<— > 



(5) 



and we can notice that by shifting the zero energy level 
in order to have (e\H\e) — 0, we can make both H ee and 
AH ee real up to the second order in time. We have thus, 
inserting eq. [5] into eq. [2l 



e**\+)+e- 



V2 



where the phases <p and v are linked by the relation: 



T 2 AH 2 eS m(v) 

4h 2 



+ o{ T 3 n 3 ) 



(9) 



The ratio between <f> (the signal we want to measure) 
and 1 — P e (the transition probability) plays for us the 
role of a signal to noise ratio (SNR), that is the quantity 
we have to maximize when we choose the operator M, 



SNR 



0(tSI) 



2(1 +cos(v)) +0{rn) 



(10) 



This quantity is not limited as a function of r and v, 
being divergent for v — > tt and r — ► with 



7T - I/ > 



(11) 



Repeating the same procedure N times on the same 
ancilla, with N of the order of !/</>, we obtain 



K , 



AK = 



t 2 AH 2 



(0) 



4^ 2 



+ 0(rW) 



where AH ee = y/ (e\H 2 \e) and f2 is a certain complex 
number whose norm is not bigger than the norm of the 
Hamiltonian's biggest matrix element divided h. 

Inserting the developments of eq. [6] into eq. [3] and IU 
we see that, due to the reality of K 2 e and AK 2 e , if we 
take v+ = ±i/_, U e ff does not alter the relative weight 
of the two eigenvectors up to the considered order. We 
will thus choose v + = — V- = v. 

The eigenvectors of M are also in this case approximate 
eigenvectors of U e f / : 



Ueff\±) 



ir 2 AH 2 am(iy) , ,=0.,, . , , 
[1 T 4 ; e 2 1 ' + 0(r 3 n 3 )}\±), (7) 



approximate, because in the error term there is a de- 
pendence upon \ip). 

The probability P e becomes 



P, 



T 2 AH 2 e (l + cos(v)) 

2h 2 



+ 0{r 6 n i ). (8) 



W{Nt+)) = \e)- 



iN<t>\ 



-iN(j>\ 



v/2 

At this point we can measure the state of the ancilla 
in the base 775G+) ± I - )), an d the probability of the 
two outcomes will be cos(N(f>) 2 and sin(N(f>) 2 . Being N(j) 
of the order of the unity, repeating the whole procedure 
a small number of times will permit us to ascertain it 
with the desired precision. Let Q be the number of mea- 
surements needed to determine this phase in the desired 
statistical confidence interval, the probability that all the 
measurements leave the system in its initial state is thus 



Nn r t 2 AHL(1 +cos(i/)) 1 ., 4 ?, h 

6 - 1 2h 2 1 

where, supposing condition [TT] to be satisfied, we can 
neglect the 0(r(l) terms . As r — > 0, 



pNQ 



Q_ 

g SNR 



that can be made arbitrarily close to one by choosing 
the right v and r. We have proved that it is indeed 
possible to measure AH ee while still impeding the system 
evolution. 
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IV. APPLICATION TO ACTUAL SYSTEMS 

With the previous procedure, we can thus measure the 
mean value of the Hamiltonian's second power in the 
initial state. Our final aim would be to determine the 
squared matrix elements that appear in the Fermi golden 
rule or, in case of reversible dynamics, the coupling fre- 
quencies with the other levels. The possibility to obtain 
these coefficients from our measurement depends on the 
details of the system and on the knowledge we have on 
the structure of the Hamiltonian. There are two simple 
cases in which we can actually succeed in calculating the 
transition frequency: first, a two levels system and, sec- 
ondly, a model in which the initial state is coupled with 
a continuum of final states and the continuum is narrow 
enough for the coupling constant not to vary significantly 
over it (in this last case, we will neglect the possible cou- 
plings between the states of the continuum, that do not 
influence AH ee ). In the first case, the Hamiltonian is 




Jt-v (10" 3 rad) 




ji-v (10" 3 rad) 



FIG. 2: The value of Pf Q , with Q = 100, is plotted using 
the exact P e (top panel) and its second order development 
(bottom panel). 



and AH ee — M7, that is we are measuring exactly 
the Rabi energy of the system. In the second case, the 
Hamiltonian is 



For r = 10~ 8 /S and it — v = 10~ 4 , the simulation gives 
a measured ratio CI/ 8 equal to one with a probability 
Pj?Q = 0.98, where the values from both panels of fig. [2] 
are = 0.99. 



H = J2 ^i\9i)(9i\ + hV\gi)(e\ + hV\e)( 9i \ (13) 

i 

and AH ee = ^|V|5^ i , that is we are measuring the 
desired squared matrix element times the integral of the 
mode density over the whole continuum. If we relax the 
hypothesis of a narrow continuum and we permit the 
coupling constant to depend on the final state, we will 
measure the quantity AH ee = ^\/Sil^l 2 ' that is the 
square root of the integral of the squared coupling con- 
stant over the continuum. From this quantity, it is not a 
priori possible (unless we know the form of the coupling 
as a function of the final state up to an overall scaling 
factor) to determine the quantity that enters the Fermi 
golden rule. This derives from the fact that, due to the 
small times implicated, the system actually sees all the 
final states, not only the resonant ones. 

The case of a two levels system (eq. H"2j) has the ad- 
vantage of allowing simple analytical and numerical cal- 
culation, thus giving us the possibility to check our pro- 
cedure, by performing calculations with the exact value 
of P e (i.e. using eq. [4] instead of eq. [8]) and to simulate 
the exact dynamics of the system and ancilla (i.e. using 
eq. instead of eq. [7]). 

In fig. H it is shown a plot of P*® for Q = S,Q = 100, 
N = l/4> with <fi calculated from eq. [5] and using re- 
spectively the exact P e defined in eq. [4] (top panel) and 
its development of eq. [5] (bottom panel). In the region 
where condition llll is verified the accord is clear. Dynam- 
ical simulations also show a good accord with the theory. 



V. CONCLUSIONS 

It has been proved that, by carefully entangling a 
metastable system with an environment, it is actually 
possible to measure a certain quantity, linked with the 
time evolution of the system, while protecting the sys- 
tem from the evolution by means of the quantum Zeno 
effect and of the entanglement itself. 

The measuring procedure here described bears some 
resemblance to the theory of protective measurements 
(@> B3 Uli) 5 at least for the fact that the wavefunction 
of the system is protected while measured, and it is thus 
possible to perform multiple measurements on the same 
quantum state. 

It can be interesting to notice that, being the quantity 
to be measured the phase of a certain effective evolution 
operator, the whole measuring procedure can be quite 
naturally restated in the language of quantum circuits by 
using a quantum phase estimation algorithm ( [H], [l3j] ) . 
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